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ABSTRACT 

Tyrosine family recombinases (YRs) are widely 
utilized in genome engineering systems because 
they can easily direct DNA rearrangement. Cre 
recombinases, one of the most commonly used 
types of YRs, catalyze site-specific recombination 
between two loxP sites without the need for 
high-energy cofactors, other accessory proteins or 
a specific DNA target sequence between the loxP 
sites. Previous structural, analytical ultracentrifuge 
and electrophoretic analyses have provided details 
of the reaction kinetics and mechanisms of Cre 
recombinase activity; whether there are reaction 
intermediates or side pathways involved has been 
left unaddressed. Using tethered particle motion 
(TPM), the Cre-mediated site-specific recombin- 
ation process has been delineated, from beginning 
to end, at the single-molecule level, including the 
formation of abortive complexes and wayward 
complexes blocking inactive nucleoprotein com- 
plexes from entering the recombination process. 
Reversibility in the strand-cleavage/-ligation pro- 
cess and the formation of a thermally stable 
Holliday junction intermediate were observed 
within the Cre-mediated site-specific recombination 
process. Rate constants for each elementary step, 
which explain the overall reaction outcomes under 
various conditions, were determined. Taking the 
findings of this study together, they demonstrate 
the potential of single-molecule methodology as 
an alternative approach for exploring reaction 
mechanisms in detail. 

INTRODUCTION 

Site-specific recombinases are categorized into the tyrosine 
family (YR) and the serine family (SR). These site-specific 
recombinases can integrate viral genomes into and excise 



them from host chromosomes. Moreover, they play a role 
in resolving multimeric chromosomes to ensure faithful 
partitioning, to maintain the correct 2\i plasmid copy 
number, and to turn a specific gene on or off (1-3). 
These recombinases possess the capability to catalyze the 
exchange of DNA strands, but their protein sequences do 
not share significant similarity (4). Based on sequence 
alignment, only five amino acid residues are conserved in 
the catalytic domains of YRs (5). It has been suggested 
that the similar catalytic properties of YRs are a result of 
their high degree of structural conservation, a similar dis- 
tribution of a-helices and P-sheets and the conservation of 
residues within the catalytic domain (6). Unlike other YR 
family recombinases that require accessory proteins and 
their cognate binding sites, Cre-/oxP and Flp-FRT have 
become useful engineering tools because of their simplicity 
(4). Both Cre-/oxP and Flp-FRT have been widely used to 
generate conditional knockouts to study gene function 
in vivo (7,8). 

YRs have been extensively studied using a variety of 
experimental approaches. For example, the recombination 
process mediated by X-Int att and Cre-/oxP has been 
studied using electron microscopy (9,10), electrophoretic 
mobility shift assay (EMSA) (11-13), atomic force micros- 
copy (14,15), analytical ultracentrifuge (AUC) (12), struc- 
tural analysis (13,16,17), surface plasmon resonance (SPR) 
(18) and single-molecule detection (19,20) analyses. 
Although Cre-/oxP is the simplest system within the YR 
family that can catalyze excisions, inversions and trans- 
locations, the reaction pathway involved is complicated. 
It has been proposed that the recombination process 
consists of four main stages. These stages are the follow- 
ing: first, two recombinase monomers sequentially bind to 
half of an inverted repeat sequence (21-23); second, two of 
these nucleoprotein complexes (presynaptic complexes) 
form one synaptic complex via a Cre-Cre interaction 
(13,24-26); third, strand cleavage and relegation are 
executed twice, which is accompanied by a conformational 
change in the recombinase molecules (4); and, finally, the 
nucleoprotein complex decomposes into separate recom- 
binases and recombinant products (27,28). Although the 
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kinetics of each step have been extensively studied using 
EMSA (29-31), AUC (12) and SPR (18), whether there 
are other side pathways and intermediates remains 
unknown. 

Recent single-molecule experiments performed by 
Mumm et al. (19) using tethered particle motion (TPM) 
have offered new insights into X Int-mediated site-specific 
recombination events from start to finish. These experi- 
ments have shown that the formation of stable X Int 
DNA complexes requires both integrase and the accessory 
proteins integration host factor and excisionase. It was 
also observed that there are side pathways and different 
nucleoprotein complexes present that are able to prevent 
inactive synaptic complexes from entering the catalytic 
recombination process. TPM has been demonstrated to 
be a powerful method that can discriminate DNA length 
topological changes with a resolution of a few base pairs 
(32,33). In this study, TPM was adapted to monitor the 
interaction between an individual DNA molecule contain- 
ing two repeated loxP sites and Cre recombinase from 
beginning to end. This detailed kinetic analysis explains 
the mechanistic events involved and their causes under 
different reaction conditions. 



MATERIALS AND METHODS 

Proteins and DNA substrates 

The phage PI site-specific recombinase Cre was purchased 
from New England Biolabs (NEB) and used without 
further purification. Two Cre mutants, Y324F and 
K201A, gifts from Professor Van Duyne, were expressed 
and purified as previously described (34). Parallel and 
inverse loxP site plasmids were constructed with pFOSl 
(NEB) using sequential restriction enzyme digestion 
(EcoRI and Xhol, NEB) according to the NEB instruc- 
tion manual. Then, 1267 bp DNA sequences containing 
either parallel loxP sites or inverse loxP sites were 
prepared via PCR using the appropriate plasmid construct 
as a template and digoxigenin and biotin-labeled primers. 
Control DNA molecules with a size of 1267 bp without 
loxP sites were prepared via PCR using pBR322 (NEB) 
as the template. DNA sequences 1267 or 580 bp in length 
containing single loxP sites were obtained via PCR using 
pFOSl as the template. The primer sequences used in 
these assays are listed in Table 1. 

Single-molecule TPM measurements and data analysis 

Streptavidin-coated polystyrene beads and the reaction 
chamber for these experiments were prepared as previ- 
ously described (32). The Brownian motion (BM) of the 
tethered beads was observed using an inverted optical 
microscope (IX-71, Olympus, lOOx objective with 
NA = 1.40) with a differential interference contrast 
(DIC) imaging system at 22°C. To enhance the contrast, 
a mercury lamp was used. Experimental images were 
acquired with an analog camera (Newvicon-70, frame 
interval: 33 ms) using custom software written in 
Labview. The centroid positions of chosen tethered 
beads were determined via a 2D Gaussian distribution. 



The subsequent 40 frames were used to calculate the x—y 
distribution scattering by means of the standard deviation 
(SD) to represent the BM amplitude. The 40-frame 
averaging window resulted in a time resolution of 1.33 s 
under our experimental setup. The data analysis followed 
previously described procedures (32). The relationship 
between the BM amplitude and the tethered DNA 
length is linear within ~1.5 kbp using a 40-frame avera- 
ging window, as shown in Supplementary Figure SI. 
The resolution in the length determination for a 1267 bp 
DNA sequence with a 40-frame averaging window is 
~150bp at a 95% confidence level (CL). A BM time 
trace was smoothed with a five-point adjacent averaging 
algorithm. Only molecules with a BM amplitude between 
70 and 91 nm with a 40-frame averaging window prior to 
Cre addition were chosen as bare 1267 bp DNA molecules 
and were then used for further analyses. Molecules that 
stuck to the glass surface (with a BM of < 26 nm) for more 
than 5 data points or that exhibited distorted movement 
(the ratio of the BM for the x coordinate to the BM for the 
y coordinate is outside the range of 0.8-1.2) at any time 
were excluded from the analyses. The dwell time distribu- 
tion histogram is fitted with either a single exponential 
or a bi-exponential decay algorithm with the following 
formulas: 

y=A x -x e ( ~ /flx " 

y = A i x e ( - k ' x 0 +A 2 x e { " k ' 2 x 0 

Ai, A 2 , k\ and /c2 are the fitting parameters determined 
by Origin 8.0. 

DNA molecules containing both parallel and inverse 
loxP sites were labeled with digoxigenin and biotin at 
their 5'-ends. The DNA molecules were immobilized on 
an anti-digoxigenin-coated coverslip; then, the other end 
of each DNA molecule was tethered to streptavidin- 
coated 200 nm polystyrene beads. Cre-mediated 
site-specific recombination was initiated by adding the 
complete reaction solution containing 50 nM Cre 
prepared in reaction buffer (50 mM Tris-HCl, pH = 7.5, 
33 mM NaCl, 10 mM MgCl 2 , 5mM dithiothreitol and 
1 mg/ml BSA). The recombination reaction was 
quenched with a 0.05% sodium dodecyl sulfate (SDS) 
solution after 30min of reaction time. To verify the 
divalent ion, pH and salt concentration dependence of 
the recombination reaction, it was performed under dif- 
ferent reaction conditions (Ca 2+ replacement, divalent ion 
deprivation, high salt and high pH). All experiments were 
performed at room temperature (22° C). 

RESULTS 

The Cre-mediated site-specific recombination targeting 
parallel loxP sites 

It has been shown that the Cre recombinase specifically 
binds to the loxP site and that formation of the synaptic 
complex includes multiple initiation steps (13,24-26). 
Once bound to the loxP sites, the Cre recombinase 
molecules bend the DNA strand sharply. Additional 
protein-protein interactions between the Cre recombinase 
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Table 1. Primer sequences used to obtain the experimental DNA sequences 



T~lN A >;i i Ktt fa tp 

l .y . \ I \ 3LI U5 11 d.LC 


Tp m til a t f 

l Gill UlclLC 


1 1 111 ill >Ll| till ill 


Parallel LoxP sites: 1267 dsDNA 


pFOSl(/rara/fe/) 


5'-DigN-GACGTGGTTTGATGGCCTC 
5'-Bio-CAGGACTGTCATTTGAGGG 


Inverse LoxP sites: 1267 dsDNA 


pFOSl(f'«V(?™») 


5'-DigN-GACGTGGTTTGATGGCCTC 
5'-Bio-CAGGACTGTCATTTGAGGG 


Single LoxP site: 1267 ds DNA 


pFOSl 


5'-DigN-GACGTGGTTTGATGGCCTC 
5'-Bio-AATAAGATCACTACCGGG 


1267 dsDNA 


pBR322 


5'-DigN-CGTCACCCTGGATGCTGTAG 
5'-Bio-CGATGGATATGTTCTGC 


580 dsDNA 


pFOSl(paralIel) 


5'-DigN-TGGTGAAAACGGGGGCG 
5'-Bio-AGGACTGTCATTTGAGGG 



molecules then stabilize the synaptic complex. Because the 
binding event between Cre recombinase and the loxP sites, 
which is followed by the inversion/excision processes, 
leads to apparent DNA length changes, a single-molecule 
TPM experiment was developed in this study to monitor 
Cre-mediated site-specific recombination in detail. 
Previous studies have demonstrated that TPM is a 
powerful tool that can be used to visualize DNA length 
and topological changes (32,33). Both continuous DNA 
length changes and conversions between different topo- 
logical states can be monitored using the nearly force-free 
TPM method (35-39). In the present study, the BM amp- 
litude of the 1267 bp DNA molecules without loxP sites 
(81.6 ± 8.6nm, n = 53, as shown in Figure 1A) was 
measured initially, which served as a control. Only a 
small proportion of molecules, 9.4%, exhibited a 
decrease in BM amplitude after a 1-h incubation with 
200 nM Cre recombinase, indicating a weak, non-specific 
interaction between Cre and the DNA (81.2 ± 9.6 nm and 
62.9 ± 11.5nm, n = 53), which is shown in Figure IB as 
well as the BM time course (Figure 1C). Next, a similar 
experiment was performed on the 1267 bp DNA molecules 
containing a single loxP site (80.5 ± 12.0 nm, n = 85, as 
shown in Figure ID). Only 15.3% of the DNA molecules 
showed a decrease in BM amplitude in response to Cre 
recombinase (80.5 ± 12.5 nm and 58.3 ± 9.4 nm, n = 85, 
Figure IE). The 1267 bp DNA molecules with a single 
loxP site were analyzed under the same scenario; the 
BM time course seldom exhibited any change during 
the 30min of reaction time (Figure IF). It is expected 
that the interaction between Cre recombinase and DNA 
should lead to a decrease in DNA length because it has 
been reported that the overall bending of the DNA strand 
in the Cre-/oxP synaptic complex occurs at an angle of 
~100° between the DNA arms (13). However, the BM 
histogram and BM time trace showed that the single inter- 
action of Cre-/oxP is too small to be discriminated by 
TPM because the resolution in the length determination 
for 1267 bp DNA molecules with a 40-frame averaging 
window is ~150bp in the range of the 95% CL 
(Supplementary Figure SI). 

Next, the 1267 bp DNA molecules containing two 
parallel loxP sites separated by 653 bp were investigated 
to explore excision recombination. Previous studies have 
suggested that two Cre recombinases bind to one of the 
loxP sites, which consists of two recombinase-binding 



elements (RBEs), and that this process occurs with a 
high level of cooperation to form the presynaptic 
complex (13,22). One such Cre-/oxP presynaptic complex 
then interacts with another Cre-/oxP presynaptic complex 
to form a recombinant synaptic complex (13,26). The BM 
value of the excision product obtained using the 580 bp 
DNA molecule was determined (43.3 ± 6.8 nm, n = 103, 
Figure 2A). Examples of single-molecule time traces are 
shown in Figure 2B(i)-(iii). To confirm whether recombin- 
ation proceeded or not, the DNA molecules were 
challenged with 0.05% SDS to disrupt non-covalent 
Cre-/oxP interactions. The change in the BM amplitude 
of the DNA molecules containing parallel loxP sites can 
be categorized as either a synapse state or an abortive 
complex state which is deficient at the Cre-Cre interaction 
interface in response to the Cre recombinase. After the 
addition of SDS, there are two possible outcomes for mol- 
ecules that have entered the synapse state. These outcomes 
are, first, a return to the original DNA substrate for the 
molecules that failed to proceed with the catalytic recom- 
bination process, as depicted in Figure 2C(i), and second, 
remaining at a short DNA length status, which corres- 
ponds to either the excision recombinant product or a 
stable Holliday junction intermediate, as depicted in 
Figure 2C(ii). In contrast, those molecules trapped at 
abortive complexes eventually returned to the original 
DNA substrate before the SDS challenge, as depicted in 
Figure 2C(iii). Among the 143 molecules analyzed with a 
starting BM amplitude between 70 and 91 nm based on 
a 40-frame averaging window, only 76 molecules showed a 
significant change in the BM time trace in response to the 
addition of Cre recombinase. Among these 76 molecules, 
15 molecules showed a decrease in the BM amplitude to 
the average value of 59.3 ± 6.8 nm within the 95% CL and 
failed to synapse within the duration of the experiment 
[Figure 2B(iii) and 2D(ii)]. This magnitude of decrease in 
the BM amplitude is different from that observed for the 
excision recombinant product obtained using 580 bp DNA 
molecules. This scenario indicates the formation of 
abortive complexes in which two loxP sites are occupied, 
but with a defect at the Cre-Cre interaction interface. 
The remainder of the molecules showed a decrease in the 
BM amplitude to an average value of 43.8 ± 8.7 nm within 
the 95% CL, similar to that of the expected excision 
recombinant product, a 580 bp DNA molecule. This ob- 
servation indicates the formation of synaptic complexes 
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Figure 1. Tethered particle motion experiment investigating the specific interaction between Cre recombinase and loxP sites. (A) Control experiment: 
the histogram of the BM amplitude was obtained using 1267 bp DNA molecules without /o.\P sites. The average BM value peaks at 81.6 ± 8.6 nm 
(n = 53). (B) A histogram of the BM amplitude obtained from 1267 bp DNA molecules without a loxP site after a 1-h incubation with 200 nM Cre 
recombinase. The histogram was fitted to a two-Gaussian distribution that represented two conditions, in which the centers of the peaks occur at BM 
values of 81.2 ± 9.6 nm and 62.9 ± 11.5nm for naked DNA and Cre-bound DNA, respectively (n = 53). (C) The change in the BM amplitude of a 
single 1267 bp DNA molecule without a loxP site in response to the addition of Cre recombinase. The dashed line indicates the addition of Cre 
recombinase. (D) Control experiment: the histogram of the BM amplitude was obtained using 1267 bp DNA molecules containing one loxV site. The 
average BM value peaks at 80.5 ± 12.0 nm (n = 85). (E) A histogram of the BM amplitude obtained from 1267 bp DNA molecules containing one 
/oxP site after 1-h incubation with 200 nM Cre recombinase. The histogram was fitted to a two-Gaussian distribution that represented two condi- 
tions, in which the centers of the peaks occur at BM values of 80.5 ± 12.5 and 58.3 ± 9.4 nm for naked DNA and Cre-bound DNA, respectively 
(n = 85). (F) The change in the BM amplitude of a single 1267 bp DNA molecule containing one loxP site in response to the addition of Cre 
recombinase. The dashed line indicates the addition of Cre recombinase. The cartoons shown below illustrate the experimental procedure. The white 
ball represents the streptavidin-labeled 200 nm polystyrene bead, the white spots represent the Cre recombinase molecules, the black curve represents 
the DNA molecules, the dashed arrowed line represents the Brownian motion amplitude, and the white line represents the loxV sequence. The error is 
the 95% CL. 



or the excision recombinant product [Figure 2D(ii)]. After 
30min of reaction time, the BM of the DNA molecules 
that failed to proceed in the catalytic recombination 
process returned to the starting BM value after the SDS 
challenge [Figure 2D(iv)]. However, the BM value of some 
DNA molecules (71.6%) remained low after the SDS chal- 
lenge, indicating the formation of either excised products 
or a stable Holliday junction intermediate [Figure 2D(iv)]. 



The dwell times between the addition of Cre and the 
shortening of the DNA molecules containing parallel 
loxP sites were classified into two groups: one for mol- 
ecules that will be trapped in an abortive complex and 
the other for molecules that will form the synaptic 
complex quickly with a short-lived presynaptic complex. 
Then, the dwell time histograms were fitted to a single ex- 
ponential model (Figure 2E and 2F) with association rate 
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Figure 2. Complex formation and recombination events of 1267 bp DNA molecules containing parallel loxP sites. (A) A histogram of the BM 
amplitude obtained from the 580 bp DNA molecules representing the recombinant products. The average BM value peaks at 43.3 ± 6.8 nm (n - 103). 
(B) (i)-( n i) The change in the BM amplitude of the 1267 bp DNA molecules containing parallel /oxP sites in response to the addition of Cre 
recombinase. (i) An example of a molecule that synapsed but failed to continue through the recombination process, (ii) An example of a molecule 
that synapsed and continued through the strand-cleavage step to the strand-migration step and then to the strand-ligation steps, (iii) An example of a 
molecule that failed to synapse within the duration of the observations. The dashed lines indicate the addition of Cre recombinase and 0.05% SDS 
after 30 min of incubation time. The bar with the punctuate pattern labels the average BM value of the expected recombinant excision product. (C). 
Reaction scheme of the Cre-mediated site-specific recombination process involving 1267 bp DNA molecules containing parallel /o.vP sites. The initial 
substrate is the 1267 bp DNA molecule containing parallel loxP sites and has an average BM amplitude of 80.5 ± 10.5 nm. Two complexes are 
formed after interaction with Cre recombinase. The first is the synaptic complex that corresponds to the formation of a Cre tetramer with an average 
BM value of 43.8 ± 8.7 nm. The other complex is the abortive complex that failed to synapse within the duration of the observations, corresponding 
to the association of Cre dimers with /oxP sites with an average BM value of 59.3 ± 6.8 nm. The reaction was stopped when 30u,l of 0.05% SDS was 
added after 30 min of incubation time. Three different outcomes were observed: (i) a molecule that synapsed but failed to complete the recombination 
process; (ii) a molecule that synapsed and formed either the recombinant product or a stable Holliday junction intermediate; and (iii) a molecule that 
failed to synapse within the duration of the observations. (D) (i) The distribution of the BM amplitude before the addition of Cre recombinase, which 
shows an average value of 80.5 ± 10.5 nm, is indicated with the letter (c). (ii) The distribution of the BM amplitude in response to the addition of Cre 
recombinase [59.3 ± 6.8 nm and 43.8 ± 8.7 nm for the abortive complex state and the synapse state, marked with (b) and (a), respectively], (iii) The 
distribution of the BM amplitude after 30 min of incubation time following the addition of Cre recombinase. (iv) The distribution of the BM 
amplitude in response to the SDS challenge after 30 min of incubation time (« = 76). The distribution of the dwell times between recombinase 
addition and the change in the BM amplitude (E) to a value representing the abortive complex state and the association rate constant of 
(6.8 ± 0.4) x 10 4 M~'s~' (R 2 = 0.96, n = 16) were obtained. (F) To a value representing the synapse state, and an association rate constant of 
(6.1 ± 0.6) x 10 4 M-' s _1 (R 2 = 0.99« = 70) was obtained. The data were fitted to a single exponential decay algorithm (Origin 8.0). The N 
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constants of (6.8 ± 0.4) x 10 4 M _1 s" 1 and (6.1 ± 0.6) x 
10 4 M _1 s _1 for the formation of an abortive complex 
and the synapse state, respectively. A total of 15 of 76 
DNA molecules containing parallel loxP sites were 
trapped in the abortive complex [Figure 2G(i)]. The dis- 
tribution of the durations of abortive complexes was fitted 
to a single exponential model, and an abortive complex 
dissociation rate constant of (5.2 ± 1.0) x 10~ 3 s _1 was 
obtained (Figure 2H). For those molecules exhibiting 
BM transitions from the original substrate to the 
synapse state [Figure 2G(ii)], leading to unsuccessful re- 
combination, the duration of the short-lived presynapse 
state is classified as type a, representing the formation of 
unstable synaptic complexes. For molecules exhibiting 
BM transitions from the synapse state to the original sub- 
strate, the duration of the short-lived presynapse state is 
classified as type b, representing the dissociation of 
short-lived presynaptic complexes. For molecules present- 
ing decreased BM, even after an SDS challenge, indicating 
the formation of stable synaptic complexes, the duration 
of the short-lived presynapse state is classified as type c, 
which represents the formation of stable synaptic 
complexes. As a result, all the durations of the short-lived 
presynapse states were categorized and used to construct 
separate histograms with a bin size of 1 s [Figure 2I(i-iii)]. 
The histograms were fitted to a single exponential model 
with rate constants of (7.7 ± 0.1) x 10" 1 s" 1 , (6.6 ± 0.4) x 
10"' s" 1 and (4.2 ± 0.7) x 10" 1 s" 1 , representing the for- 
mation of unstable synaptic complexes, the decomposition 
of unstable presynaptic complexes and the formation 
of stable synaptic complexes of DNA molecules that 
successfully synapse during the experimental period, 
respectively. 

The Cre-mediated site-specific recombination targeting 
inverse loxP sites 

To distinguish whether Holliday junction intermediates or 
excision products are formed within the 30-min incubation 
time, the same experiments were performed on 1267 bp 
DNA molecules containing inverse loxP sites separated 
by 653 bp. The recombination process in the DNA mol- 
ecules containing inverse loxP sites results in an inverted 
653 bp DNA segment between the two sites and an un- 
changed DNA length for the recombinant product. 



Examples of single-molecule time traces are shown in 
Figure 3A(i)-(iii). For molecules that have proceeded to 
the synapse state, there are two possible outcomes fol- 
lowing an SDS challenge: first, a return to the original 
DNA substrate or the formation of an inverse recom- 
binant product [Figure 3B(i)], and second, remaining as 
a shortened DNA product, representing the stable 
Holliday junction intermediate [Figure 3B(ii)]. Molecules 
that are trapped at the abortive complex will return to the 
original DNA substrate prior to the SDS challenge. Of the 
301 analyzed molecules with a starting BM value of 
79.8 ± 7.5 nm, 101 molecules exhibited a change in their 
BM amplitude in response to the addition of Cre recom- 
binase. Of the 101 analyzed molecules, 18 showed a 
decrease in BM amplitude to an average value of 
64.3 ± 3.6 nm within the 95% CL and failed to synapse 
within the experimental duration [Figure 3A(iii) and 
3C(ii)], indicating the formation of abortive complexes. 
When the remaining analyzed molecules were examined, 
the BM amplitude was found to decrease rapidly to 
an average value of 48.6 ± 10.2 nm within the 95% 
CL, indicating the formation of synaptic complexes 
[Figure 3C(ii)]. The BM amplitude of the synapse state 
for the DNA molecules containing inverse loxP sites is 
different from that of DNA molecules containing 
parallel loxP sites. The discrepancy in the BM amplitude 
between the synapse states occurs because the two loxP 
sites, when in the correct synapse state, are expected to be 
in an anti-parallel orientation. This constraint results in 
extra topological bending of DNA molecules containing 
inverse loxP sites (13,16,40). After 30min of incubation 
with Cre recombinase, the percentage of DNA molecules 
in the abortive complex decreased because they were 
returning to the original DNA substrate state 
[Figure 3C(iii)]. Immediately upon the SDS challenge at 
30min reaction time, the BM amplitude of most of the 
DNA molecules (90%) remained at a low value of 
48.6 ± 10.2 nm, indicating the formation of stable 
Holliday junction intermediates. Similar association con- 
stants of (8.1 ± 0.5) x 10 4 M-'s-' and (6.4 ± 0.3) x 
10 4 M _1 s _1 for the abortive complexes and synaptic 
complexes, respectively, were calculated by fitting 
the dwell histograms to a single exponential model 
(Figure 3D and E). 



Figure 2. Continued 

mentioned above is the number of molecules that were observed. The error is within the 95% CL. (G) BM time-trace of a 1267 bp DNA molecule 
containing parallel loxP sites in response to the addition of Cre recombinase. The dashed lines indicate the addition of Cre recombinase and the 
addition of 0.05% SDS after 30min of incubation time. The bar with the punctuate pattern indicates the average BM value of the expected excision 
recombinant product, (i) The BM time trace of a molecule that has successfully formed a synaptic complex following a short-lived presynapse state, 
(ii) The BM time trace of a molecule that failed to synapse and was trapped in the abortive complex state. For those molecules undergoing BM 
transitions from the original substrate to the synapse state, leading to an unsuccessful recombination, the short-lived presynapse state duration is 
classified as type a. For molecules that experience BM transitions from the synapse state to the original substrate, the short-lived presynapse state 
duration is classified as type b. For those molecules showing decreased BM, even after the SDS challenge, indicating the formation of a stable 
synapse state, the short-lived presynapse state duration is classified as type c. (H) The dwell times in the abortive complexes formed with the DNA 
molecules containing parallel loxP sites that failed to synapse during observation were pooled and then fitted to a single exponential decay 
[R 2 = 0.98, rl = (5.2 ± 1.0) x 10~ 3 s~', n = 26]. (I) (i)-(iii) The dwell times in the short-lived presynapse state were pooled separately to build 
dwell-time histograms fitted to a single exponential decay algorithm with values of (7.7 ± 0.1) x 10~ s~ (n = 17, R 2 = 0.99), (6.6 ± 0.4) x 10~' s -1 
(n = 16, R 2 = 0.97) and (4.2 ± 0.7) x 10~' s (;; = 40, R 2 = 0.99), representing the association rate constant of unstable synaptic complexes, the 
dissociation rate constant of presynaptic complexes and the association rate constant of correct synaptic complexes, respectively. The N mentioned 
above is the number of events observed. The error is within a 95% CL. 
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Figure 3. Complex formation and recombination events of 1267 bp DNA molecules containing inverse loxP sites. (A) (i)-(iii) The change in the BM 
amplitude of the 1267 bp DNA molecules containing inverse loxP sites in response to the addition of Cre recombinase. (i) An example of a molecule 
that synapsed and then returned to the original substrate, corresponding to either no reaction or completion of recombination, (ii) An example of a 
molecule that synapsed and was then trapped in a stable Holliday junction intermediate, (iii) An example of a molecule that failed to synapse within 
the duration of the observations. The dashed lines indicate the addition of Cre recombinase and 0.05% SDS after 30min of incubation time. The bar 
with the punctuate pattern indicates the average BM value of the expected excision recombinant product. (B) Reaction scheme of the Cre-mediated 
site-specific recombination process involving 1267 bp DNA molecules containing inverse loxP sites. The starting substrates are 1267 bp DNA mol- 
ecules containing inverse loxP sites and exhibit an average BM amplitude of 79.8 ± 7.5 nm. Two complexes are formed after interaction with Cre 
recombinase. The first complex is the synaptic complex, which corresponds to the formation of a Cre tetramer with an average BM value of 
48.6 ± 10.2 nm. The other complex is the abortive complex that failed to synapse within the duration of the observations and corresponds to the 
association of Cre dimers with loxP sites with an average BM value of 64.3 ± 3.6 nm. The reaction was stopped when 30 ul of 0.05% SDS was added 
after 30min of incubation time. Three different outcomes were observed: (i) a molecule that synapsed but either failed to complete recombination or 
formed the recombinant product after completion of the recombination process; (ii) a molecule that synapsed but was trapped in a stable Holliday 
junction intermediate; and (iii) a molecule that failed to synapse within the duration of the observations. (C) (i) The distribution of the BM amplitude 
before the addition of Cre recombinase showed an average value of 79.8 ± 7.5 nm, which is indicated with (c). (ii) The distribution of the BM 
amplitude in response to the addition of Cre recombinase (64.3 ± 3.6, 48.6 ± 10.2 nm for the abortive complex state and the synapse state, marked 
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Kinetic analysis and the detailed reaction mechanism of 
Cre-mediated site-specific recombination 

For the 1267 bp DNA molecules containing parallel loxP 
sites, 25 of 76 molecules showed fluctuating BM transi- 
tions between the synapse state and the original substrate 
within the observation period (Figure 4A). Because 
excised DNA molecules are the expected recombinant 
products of DNA molecules containing parallel loxP 
sites, the synapse state duration is correlated with the sta- 
bility of unstable synaptic complexes, i.e. wayward 
complexes. Therefore, the dwell time distribution can be 
interpreted in terms of single exponential decay. Similar 
fluctuation phenomena were observed for DNA molecules 
containing inverse loxP sites; the synapse state duration 
reflects both the stability of the wayward complex and 
the catalytic reactivity of recombination. Thus, in these 
circumstances, the synapse state duration distribution 
can be interpreted in terms of a bi-exponential model. 
Surprisingly, the duration histograms for both the 
parallel and inverse loxP site-containing DNA molecules 
could be fitted to a bi-exponential model, and similar rate 
constants were obtained (Figure 4B and Supplementary 
Figure S2 fitted via least square and maximum likelihood 
estimation, respectively). For the Cre-parallel loxP sites 
DNA system, the bi-exponential decay behavior can 
only be observed when the recombination process, i.e. 
strand cleavage/ligation, is reversible. In our experiments, 
the Holliday junction intermediate was very stable and 
served as a reaction trap, and the reaction could go 
forward, to form excision recombinant products, or go 
backward to the original substrates with an equal prob- 
ability. The fast reaction process, with rate constants of 
(2.7 ± 0.5) x 10" 2 s _1 and (2.7 ± 0.6) x 10" 2 s _1 for DNA 
molecules containing parallel and inverse loxP sites, re- 
spectively, represents the dissociation rate constant of 
the wayward complex. The slow reaction process, with 
rate constants (2.0 ± 0.1) x 10" 3 s" 1 and (2.2 ± 0.2) x 
10~ 3 s _1 for the DNA molecules containing parallel and 
inverse loxP sites, respectively, represents the catalytic re- 
activity of recombination. 

To verify whether the synapse state duration for the 
DNA molecules containing two loxP sites reflects both 
the decomposition of the wayward complex and the com- 
pletion of recombination and whether cleavage is required 
for the formation of stable synaptic complexes, a Cre 
mutant (Y324F) was used to perform the same experi- 
ments and analyses. As expected, Tyr 324 is necessary 
for recombination because it serves as a nucleophile that 
attacks the phosphodiester bond within the active site of 
the synaptic complexes. Based on a previous study, it was 
noted that none of the 19 replacement residues showed 
any detectable recombination activity (41). Moreover, it 



has been shown that subsequent DNA cleavage by Int, 
a membrane of the tyrosine site-specific recombinases, 
will stabilize the synaptic complex; otherwise, the 
synaptic complex will easily decompose (19). Typical 
BM time traces are shown in Figure 5A(i-ii). The 
anticipated reaction scheme is depicted in Figure 5B, 
where the change in BM amplitude in response to 
CreY324F can be classified into either a synaptic 
complex or an abortive complex. As expected, the defi- 
ciency in strand cleavage will lead to no recombinant 
product being formed, as verified by the SDS challenge 
[Figure 5B(i)], and those molecules that formed abortive 
complexes will eventually return to the original DNA sub- 
strate [Figure 5B(ii)]. Among the 85 molecules with a 
starting BM amplitude between 70 and 91 nm, 42 mol- 
ecules exhibited a significant change in the BM amplitude 
in response to the addition of CreY324F recombinase. A 
total of 7 of 42 molecules exhibited a decrease in the BM 
amplitude to an average value of 62.0 ± 4.2 nm and failed 
to synapse within the experimental period [95% CL, 
Figure 5C(ii)], indicating the formation of abortive 
complexes. For the remaining molecules, the BM ampli- 
tude decreased to an average value of 47.5 ± 9.4 nm, and a 
short-lived presynapse state was observed [95% CL, 
Figure 5C(ii)], indicating the formation of synaptic 
complexes. After 30min of reaction time, most of the 
molecules returned to the original DNA substrates 
[Figure 5C(iii)]. After the SDS challenge, all the molecules 
were expected to return to the original DNA substrates 
because CreY324F lacks the ability to carry out strand 
cleavage. The fact that there were two molecules trapped 
at the lower BM value could be an experimental artifact or 
indicate purification contamination or sticking problems. 
With exception of the lower percentage of experimental 
artifacts, the TPM results shown here are consistent with 
those obtained by Gibb et al. (41). Smaller association rate 
constants of (1.6 ± 0.2) x 10 4 M _1 s" 1 and (1.1 ± 0.0) x 
10 4 M~'s -1 for the formation of abortive complexes 
and synaptic complexes, respectively, were obtained by 
fitting the duration histograms to a single exponential 
model [Figure 5D(i-ii)]. An abortive complex dissociation 
rate constant of (1.1 ±0.1)x 10~ 2 s _1 was obtained by 
fitting the abortive complex duration distribution to a 
single exponential model [Figure 5E(i)]. When the remain- 
ing molecules that formed a synaptic complex rapidly, 
with a short-lived presynapse state, were analyzed 
further, all of them (excluding experimental artifacts) ex- 
hibited fluctuating transitions between the synapse state 
and the original DNA substrate [Figure 5B(ii)]. The 
synapse duration histogram was fitted to a single exponen- 
tial model with a rate constant of (3.2 ± 0.3) x 10~ 2 s _1 , 
which represented the dissociation of wayward complexes. 



Figure 3. Continued 

with (b) and (a), respectively), (iii) The distribution of the BM amplitude after 30min of incubation time following the addition of Cre recombinase. 
(iv) The distribution of the BM amplitude in response to the SDS challenge after 30min of incubation time (n — 101). The distribution of the dwell 
times between recombinase addition and the change in the BM amplitude (D) to a value representing the abortive complex state and an association 
rate constant of (8.1 ± 0.5) x 10 4 M~' s _1 (R 2 = 0.98, n = 18) were obtained. (E) To a value that represents the synapse state, and an association rate 
constant of (6.4 ± 0.3) x 10 4 M~' s" 1 (R 2 = 0.98, n = 81) was obtained. The data were fitted to a single exponential decay algorithm (Origin 8.0). The 
N mentioned above is the number of molecules that were observed. The error is within the 95% CL. 
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Figure 4. Kinetic analysis of the dwell time in the synapse state. (A) 
The BM time-trace of a 1267 bp DNA molecule containing parallel 
loxF sites in response to the addition of Cre recombinase. The 
dashed lines indicate the addition of Cre recombinase and 0.05% 
SDS after 30min of incubation time. The bar with the punctuate 
pattern indicates the average BM value of the expected excision recom- 
binant product. (B) (i) The distribution of the dwell time in the synapse 
state, indicated by the double-arrowed line in the above time trace, for 
DNA molecules containing parallel loxP sites was pooled to construct a 
histogram, and the rate constants were obtained by fitting to either a 
single exponential decay algorithm [R 2 = 0.91, k x = (4.4 ± 0.5) x 
10~ 3 s~\ n = 118] or a bi-exponential decay algorithm [ii 2 = 1 .00, 
fcj = (2.7 ± 0.5) x 10~V\ ,4, =0.80; k 2 = (2.0 ± 0.1) x 10" 3 s _1 , 
^2 = 0.20, n = 118]. Solid and dashed lines are used to represent 
single and bi-exponential fitting curves, respectively, (ii) The distribu- 
tion of the dwell time in the synapse state, indicated by the 
double-arrowed line in the above time trace, for DNA molecules con- 
taining inverse loxP sites was pooled to construct a histogram, and the 
rate constants were obtained by fitting to either a single exponential 
decay algorithm [R 2 = 0.91, k x = (6.5 ± 0.9) x 10 _3 s _1 , n = 178] or a 
bi-exponential decay algorithm [R 2 = 1.00, k x = (2.7 ± 0.6) x 10~ 2 s~\ 
A, = 0.85; k 2 = (2.2 ± 0.2) x 10~ 3 s _1 , A 2 = 0.15, n = 178]. Solid and 
dashed lines are used to represent single and bi-exponential fitting 
curves, respectively. The N mentioned above is the number of events 
observed. The error is within the 95% CL. 



The kinetic reaction scheme of CreY324F-mediated re- 
combination for DNA molecules containing parallel 
loxP sites is shown in Figure 5F. It can be inferred that 
CreY324F is deficient in the formation of a stable synaptic 
complex or lacks cleavage ability; therefore, the catalysis 



recombination pathway is blocked in the CreY324F-/oxP 
system. 

It has been known that the cleavage-deficient mutant 
CreK201A, from which the required general acid 
catalyst that activates 05' has been removed, possesses a 
similar synapsis ability as wild-type Cre, with K d = 8.9 ± 
0.6 nM (12). Therefore, CreK201A was used to perform 
the same experiments and analyses in order to verify 
whether cleavage is required for the formation of the 
stable synaptic complex in the Cre-/oxP system. The BM 
time traces are shown in Figure 6A, and there are two 
typical reaction behaviors, similar to those of CreY324F. 
One reaction behavior is that a DNA molecule formed a 
synaptic complex, but there was no formation of a recom- 
binant product verified by the SDS challenge, and the 
other is that a DNA molecule failed to synapse during 
the duration of the observations. Analysis of 59 molecules 
with a BM amplitude between 70 and 91 nm showed that 
43 molecules exhibited a significant change in BM in 
response to the addition of CreK201A recombinase. The 
BM histograms are shown in Figure 6B, which are similar 
to those observed for CreY324F. Association rate con- 
stants of (2.0 ± 0.1) x 10 4 M _1 s"' and (1.5 ± 0.1) x 
10 4 M _1 s _1 for the formation of abortive complexes 
and synaptic complexes, respectively, were obtained by 
fitting the duration histograms to a single exponential 
model [Figure 6C(i-ii)]. The stability of the abortive 
complex was obtained by fitting the abortive complex 
duration histogram to a single exponential model with 
a dissociation rate constant of (4.0 ± 0.7) x 10~ 3 s -1 
[Figure 6D(i)]. When the duration of the synaptic 
complex was analyzed further, it was found that the 
duration distribution was similar to that observed in the 
wild-type Cre-/oxP system, and a bi-exponential decay 
model was used. The fast process, with a rate constant 
of (4.1 ± 1.9) x 10 _2 s _1 , represents the decomposition of 
the wayward complex. The same rate constants were 
obtained for both wild-type Cre and CreY324F, which 
reinforced the existence of a wayward complex, in the 
Cre-/oxP system. On the other hand, it has been shown 
that the catalysis recombination process, i.e. strand 
cleavage/ligation, is reversible. Limited by the resolution 
of the TPM, the change in BM amplitude can only signal 
the decomposition of the synaptic complex, and the BM 
amplitudes observed for synaptic complexes, cleavage 
synaptic complexes and Holliday junction intermediates 
showed insignificant differences. Moreover, it has been 
reported that the rate-limiting step in the Cre-mediated 
recombination process is the dissociation of the 
recombined synapse (31). Therefore, the slow process 
observed in wild-type Cre and CreK201A, with a rate 
constant of (2.6 ± 0.2) x 10 _3 s _1 , represents the decom- 
position of the stable synaptic complex as well as the de- 
composition of the recombined synaptic complex. 

Effects of pH, divalent ions and monovalent ions on the 
Cre-mediated site-specific recombination process 

Previous studies have shown that the Cre-/oxP synapsis is 
pH dependent and that the presence of divalent ions or 
polyamines will increase the recombination efficiency by 



Nucleic Acids Research, 2012, Vol. 40, No. 13 621 7 




SDS challenge 

245bp 
J 301bp( 



B 



Synaptic complex 



BM=79.4±7.7nm 





•* 





0 



No reaction (i) 



No reaction (ii) 



Abortive complex 



Time(s) 

(a) (b) 



20 
10 
0 

20 

>> 

o 10 
c 

CD 0 

| 20 
LL 10 

0 

20 
10 
0 



0 



D 



80 



100 

(ii) 



(I 


20 


40 

. F=E 


6$ 

I i - 


i 


* 


100 

(iii) 

1 — 1 


(I 


20 


40 


60 


80 


100 














(iv) 










n 




|— Ii— , 


0 


20 


40 


60 


80 


100 




500 1000 1500 2000 

Abortive complex association time (s) 



2000 




E 

40 
>■ 30 
S 20 

3 

Q) lu 

t ° n 
> 200 Q 



Brownian motion amplitude (nm) 



(i) 



500 1000 1500 

Synapse association time (s) 



150 



5 100 
o 

< 50 



100 200 300 400 500 600 700 Presynaptic 

Abortive complex duration (s) 



, , , | A (7.1±0.5)xl0" 1 s" 1 
<l.l±0.0 4 )xlO<M^ (1.1±0.1)X10-S>^ 



Parallel loxP substrate 

(1.6+0.2JX10 4 M" 1 s 1 




(ii) 



0 100 200 300 400 500 600 700 

Synapse state duration (s) 



complex 

(5.9±0.4)x^4t( 3 - 2±0 - 3 > Xl °- 2 - 



Wayward 

complex 

P=83.3% 




Abortive 
complex 
P=16.7% 



Figure 5. Complex formation and recombination events of 1267 bp DNA molecules containing parallel loxP sites mediated by cleavage-deficient 
CreY324F. (A) (i)-(ii) A change in the BM amplitude of 1267 bp DNA molecules containing parallel loxP sites in response to the addition of 
CreY324F. (i) An example of a molecule that synapsed but failed to continue to the recombination process, (ii) An example of a molecule that 
failed to synapse within the duration of the observations. (B) The reaction scheme of the Cre-mediated site-specific recombination process in a 1267 bp 
DNA molecule containing parallel /o.\P sites. The initial substrate is a 1267 bp DNA molecule containing parallel loxP sites with an average BM value 
of 79.4 ± 7.7 nm. Two complexes are formed after the addition of Cre recombinase. The first complex is the synaptic complex, corresponding to the 
formation of a Cre tetramer with an average BM value of 47.5 ± 9.4 nm. The other complex is the abortive complex that failed to synapse within the 
duration of observations, corresponding to the association of Cre dimers with loxP sites, with an average BM value of 62.0 ± 4.2 nm. When 30 ul of 
0.05% SDS was added after 30min of incubation time, the reaction was stopped. Two different outcomes were observed: (i) a molecule that synapsed 
but lacked the strand-cleavage capability to complete recombination and returned to the original DNA substrate; and (ii) a molecule that failed to 
synapse within the duration of observations and eventually returned to the original DNA substrate. (C) (i) The distribution of the BM amplitude 
before the addition of Cre recombinase, with an average value of 79.4 ±7.7 nm, which is indicated with (c). (ii) The distribution of the BM amplitude 
in response to the addition of Cre recombinase (62.0 ± 4.2 and 47.5 ± 9.4 nm for the abortive complex state and the synapse state, marked with (b) and 
(a), respectively), (iii) The distribution of the BM amplitude after 30 min of incubation time, (iv) The distribution of the BM amplitude in response to 
the SDS challenge after 30 min of incubation time (n = 42). The distribution of the dwell times between recombinase addition and the change in BM 
amplitude (D) to a value representing the abortive complex state and an association rate constant of (1.6 ± 0.2) x 10 4 M~' s (R 2 = 0.96, n = 18) were 
obtained. (F) To a value that represents the synapse state, and an association rate constant of (1.1 ± 0.0) x 10 4 M~'s~' (R 2 = 0.99, n = 64) was 
obtained. The data were fitted to a single exponential decay algorithm (Origin 8.0). The N mentioned above is the number of molecules that were 
observed. (E) (i) The dwell times in the abortive complex state for DNA molecules containing parallel loxP sites that failed to synapse during the 
observations were pooled and then fitted to a single exponential decay algorithm [R 2 = 0.98, rl = (1.1 ± 0.1) x 10~ 2 s~', n = 40]. (ii) The distribution 



of the dwell times in the synapse state was fitted to either a single exponential decay algorithm [R 



0.98, k. 



(3.2 ± 0.3) x 10 -s , n = 189] or a 

bi-exponential decay algorithm [R 1 = 1.00, k, = (6.3 ± 0.5) x 10~ 2 s~', A, = 0.95; k 2 = (5.7 ± 0.5) x 10~ J s~', A 2 = 0.05, n = 189]. The solid and 
dashed lines represent single and bi-exponential fitting curves, respectively. The N mentioned above is the number of events that were observed. 
The error is within the 95% CL. (F) The mechanism through which bacteriophage PI CreY324F mediates site-specific recombination. The arrows on 
the DNA molecule indicate the parallel loxV sites. The empty circles represent Cre recombinases. 
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Figure 6. Complex formation and recombination events for 1267 bp DNA molecules containing parallel loxP sites mediated by cleavage-deficient 
CreK201A. (A) (i)-(ii) A change in the BM amplitude of 1267 bp DNA molecules containing parallel loxP sites in response to the addition of 
CreK201A. (i) An example of a molecule that synapsed but failed to continue to the recombination process, (ii) An example of a molecule that failed 
to synapse within the duration of the observations. (B) (i) The distribution of the BM amplitude before the addition of Cre recombinase, with an 
average value of 78.8 ± 10.2 nm, which is indicated with (c). (ii) The distribution of the BM amplitude in response to the addition of Cre recombinase 
(61.0 ± 1.3 nm and 44.0 ± 12.0 nm for the abortive complex state and the synapse state, marked with (b) and (a), respectively), (iii) The distribution 
of the BM amplitude at 30min of incubation time, (iv) The distribution of the BM amplitude in response to the SDS challenge at 30min of 
incubation time (n = 43). (C) The distribution of the dwell times between recombinase addition and the change in the BM amplitude to (i) a value 



that represented the abortive complex state and an association rate constant of (2.0 ± 0.1) x 10 M 



(ir = 0.98, n = 16) and to (ii) a value that 



represented the synapse state and an association rate constant of (1.5 ± 0.1) x 10 M 



(R = 0.98, = 121) were obtained. The data were fitted 



to a single exponential decay algorithm (Origin 8.0). The N mentioned above is the number of molecules observed. (D) (i) The dwell times in the 
abortive complex state for the DNA molecules containing parallel loxP sites that failed to synapse during the observations were pooled and then 
fitted to a single exponential decay algorithm [R 1 = 0.98, t, = (4.0 ± 0.7) x 10~ 



state was fitted to either a single exponential decay algorithm [R 1 = 0.98, k\ = (6.8 ± 0.8) x 10 s 1 



, n = 7]. (ii) The distribution of the dwell times in the synapse 
n = 148] or a bi-exponential decay algorithm 



[R 1 = 1.00, fcj = (4.1 ± 1.9) x 10~~s , A t = 0.90; k 2 = (2.6 ± 0.2) x 10~ J s , A 2 = 0.10, n = 148]. Solid and dashed lines represent single and 
bi-exponential fitting curves, respectively. The N mentioned above is the number of events observed. The error is within the 95% CL. (E) The 
mechanism through which bacteriophage PI CreK201A mediates site-specific recombination. The arrows on the DNA molecule indicate the parallel 
/oxP sites. The empty circles represent Cre recombinases. The arrows on the DNA molecule indicate the parallel /oxP sites. 



facilitating cooperative Cre-/o.\-P binding (16,18,25). 
However, the causes of the occurrence of unstable 
synaptic complex breakdown at a high pH and how 
divalent ions facilitate the recombination process remain 
unclear (12). In this study, the potential of TPM to delin- 
eate the detailed kinetic steps involved in this process has 
been demonstrated, which allows a more quantitative 
understanding of the phenomena involved. To expand 
our knowledge regarding the effects of pH and ions, 
Cre-mediated site-specific recombination was studied 
using TPM under five different reaction conditions: in 
the presence of Mg 2+ , in the presence of Ca 2+ , in the 
absence of divalent ions, in the presence of a high concen- 
tration of monovalent ions and at a high pH. The Mg 2+ 
condition was defined as the normal condition and has 
been explored above. The same experiments and 



analyses were then performed under the various condi- 
tions described above, and the resulting kinetic parameters 
are listed in Table 2. The recombination efficiency in these 
circumstances is defined as the ratio of the number of 
molecules forming either recombinant products or 
Holliday junction intermediates (with a BM amplitude 
of 43.8 ± 8.7 nm after the SDS challenge at 30min of in- 
cubation time) to the number of molecules responding to 
Cre recombinase. The detailed kinetic analyses are pre- 
sented in Supplementary Figures S3-S6. 

It is known that the catalytic activity, but not the DNA 
binding affinity of some endonucleases is affected when 
Mg 2+ ions are replaced by Ca 2+ ions (42,43). According 
to previous studies and gel electrophoresis results, it has 
been shown that there is no significant difference in the 
Cre-mediated recombination process under divalent ion 
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Figure 7. Mechanism by which bacteriophage PI Cre recombinase mediates site-specific recombination. The rate constant for each step was obtained 
via either single exponential decay fitting or bi-exponential decay fitting to the duration distribution histograms. The arrows on the DNA molecule 
indicate the parallel loxP sites. The empty circles represent the Cre recombinases. The small black dots represent the covalent phosphotyrosine 
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replacement conditions (12). Therefore, TPM experiments 
and kinetic analyses were performed using different reac- 
tion buffer conditions to explore the effect of two divalent 
ion species, Mg 2+ and Ca 2+ , on the Cre-mediated recom- 
bination process; in addition, the effect of the absence of 
divalent ions was also explored. When the TPM results 
were examined in depth, the affinity of Cre-/oxP and the 
stability of the synaptic complex were found to be similar 
in the presence of either Mg 2+ or Ca 2+ (Table 2, 
Supplementary Figure S3). These results reinforced the 
hypothesis that Mg 2+ and Ca 2+ have the same positive 
effect on the Cre-mediated recombination process (25). 
When divalent ions were removed from the reaction 
buffer, it was found that the recombination efficiency 
was lower than under normal conditions. A detailed 
kinetic analysis showed that the lower recombination effi- 
ciency could be attributed to a lower probability (29.4%) 
of forming the correct synaptic complex leading into a 
successful catalytic recombination process. In addition, 
based on kinetic analysis, it was found that not only was 
the Cre-/o.rP interaction weaker, but the stability of the 
synaptic complex was also reduced in the absence of 
divalent ions (Table 2, Supplementary Figure S4). 
Structural data indicated that the enhancement of the 



stability of the synaptic complex was due to the formation 
of divalent-cationic species stabilized by close contacts 
(~4A) between the phosphate backbones within the 
Cre-/oxP synaptic complex (12,16). Based on the experi- 
mental results obtained under divalent ion replacement 
and divalent ion deprivation, the requirement for 
divalent ion species during the Cre-mediated recombin- 
ation process has been confirmed. 

To explore the effect of high concentrations of mono- 
valent ions, the same experiments and kinetic analyses 
were performed at a concentration of 150mM Na + . The 
electrostatic interaction between the recombinases and the 
loxP site appears to be reduced at a high ionic strength; 
hence, the association rate constants have been found to 
be smaller (44,45). However, the kinetic analysis showed 
that the high concentration of monovalent ions appeared 
to only affect the formation and stability of abortive 
complexes and not the formation and stability of 
synaptic complexes (Table 2, Supplementary Figure S5). 
These results are consistent with previous studies in which 
it was found that divalent cationic species alone are suffi- 
cient to stabilize the close contacts between the phosphate 
backbones within the synaptic complexes (12). Based on 
the TPM kinetic analysis, it can be concluded that the 
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Table 2. Kinetic and thermodynamic information related to Cre-mediated site-specific recombination under different conditions 



Conditions 


[NaCl] 
mM 


[M 2+ ] 
mM 


PH 


^Abortive 

(%) 


P Wayward 

(%) 


^Association 

(M-'s -1 )" 


Normal 
(« = 76) 


33 


2 


7.5 


19.7 


33.7 


(6.8 ± 0.4) x 10 4 
(6.1 ± 0.6) x 10 4 


Ca 2+ (« = 37) 


33 


2 


7.5 


IX. 2 


31.8 


(5.9 ± 0.7) x 10 4 
(5.0 ± 0.3) x 10 4 


No divalent 
cations 
(« = 61) 


33 


0 


7.5 


17.3 


53.3 


(6.4 ± 0.5) x 10 4 
(4.4 ± 0.3) x 10 4 


High salt 

(« = 47) 


150 


2 


7.5 


13.5 


54.1 


(5.9 ± 0.7) x 10 4 
(5.0 ± 0.3) x 10 4 


pH = 8.5 
(» = 82) 


33 


2 


8.5 


17.2 


44.8 


(5.7 ± 0.1) x 10 4 
(4.0 ± 0.2) x 10 4 



^Dissociation 



n synapse 
decomposition 



^Recombined synapse Efficiency 
decomposition (,S ) ( /o) 



(2.0 ± 0.1) x 10" 3 71.6 



" short-lived 

(s-') b 

ptvsynapse V a / 

(7.7 ± 0.1) x 10~' (2.7 ± 0.5) x 10~ 2 
(6.6 ± 0.4) x 10~' 
(4.2 ± 0.7) x 10 1 

(1.9 ± 0.3) X 10 -1 (2.6 ± 0.4) x 10" 2 (2.0 ± 0.2) x 10" 3 65.7 
(2.8 ± 0.2) x 10" 1 
(6.7 ± 0.9) x 10"' 

(3.6 ± 0.3) x 10 _1 (2.2 ± 0.2) x 10~ 2 (1.7 ± 0.1) x 10" 3 
(2.3 ± 0.2) x 10~' 
t a q ^m\ v in-' 



5.4 ± 0.7) x 10 -3 
3.0 ± 0.8) x 10~ 3 



)x 10" 



(3.6 ± 0.3) x 10"' 
(2.3 ± 0.2) x 10"' 
(4.8 ± 0.2) x 10" 1 

(1.9 ± 0.3) x 10"' (1.8 ±0. 
(2.8 ± 0.2) x 10" 1 
(6.7 ± 0.9) x 10~' 

(9.6 ± 0.9) x 10~' (3.4 ± 0.9) x 10 
(9.0 ± 0.8) x 10~' 
(1.1 ± 0.1) 



.2) x 10~ 2 (2.2 ± 0.2) x 10~ 3 
(2.5 ± 0.1) x 10~ 3 



37.1 



44. S 



54.6 



^Association 

b 



: the two listed constants represent the association rate constants of the presynaptic complex and the synaptic complex. 



^short-lived P resynapse : the three listed constants represent the rate constants for the formation of unstable synaptic complexes, the decomposition of 
unstable presynaptic complexes and the formation of stable synaptic complexes of DNA molecules that successfully synapse during the experimental 
period. 



lower recombination efficiency observed at high monova- 
lent ionic strengths is due to a lower probability (32.4%) 
of forming the correct synaptic complex, rather than a 
reduction of the catalytic ability (Table 2). 

Previous studies have indicated that the Cre-mediated 
site-specific recombination reaction is robust over a pH 
range of 5.5-8.0 but is inhibited once the pH is raised to 
8.5. A possible explanation for the pH dependence of 
synapsis and recombination is that the cysteine residues 
located within the helical core binding domain of Cre re- 
combinase tend to ionize when the pH is greater than 8.33 
and that these ionized residues affect the stability of 
synaptic complexes (12). The TPM results obtained at a 
high pH showed that a lower association rate occurred, 
which appears to confirm that increased electrostatic re- 
pulsion is the main factor that weakens the Cre-/oxP inter- 
action. The kinetic analysis showed that the increased 
electrostatic repulsion at a high pH weakens the Cre- 
Cre interaction [fc Association = (4.0±0.2)xl0 4 M _1 s _1 ], 
reducing the probability of forming the correct synaptic 
complex (38%, compared to 46.6% under normal condi- 
tions) (Table 2, Supplementary Figure S5). The results 
regarding the pH dependence of synapsis obtained using 
TPM together with detailed kinetic analysis reinforced 
the hypothesis that increased electrostatic repulsion is 
the main factor that affects synaptic complex stability at 
a high pH. 

DISCUSSION 

Using the above kinetic results, the kinetic details of the 
reaction pathway from beginning to end can be delineated 
(Figure 7). In a previous single-molecule study on the 
tyrosine family recombinase X Int, it was found that pre- 
synaptic complexes formed at a rate constant >0.05s -1 , 



and no detectable dissociation process can be observed 
within 60min of reaction time (19). Moreover, it 
has been found that the interaction between Cre and 
loxP occurs very rapidly, with a rate constant of 
~10 8 M~'s _1 (31). However, the association rate 
constant of (6.8 ± 0.4) x 10 4 M _1 s" 1 obtained via TPM 
is slower compared to previous results. It is believed that 
Cre recombinase is monomeric in solution and binds to 
the loxP site one by one to form a dimer (21-23). 
However, only the presynapse state in which two loxP 
sites are occupied by Cre recombinases can be 
differentiated because of the limited resolution of TPM. 
According to our data, the ^Association value observed is the 
rate constant of abortive complex formation, rather than 
that of DNA binding. Therefore, it cannot be ruled out 
that the measured association rate constant is 
underestimated. For these abortive complexes, the dissoci- 
ation rate constant of Cre recombinase from the loxP site 
is (5.2 ± 1.0) x 10~ 3 s _1 , which is consistent with previous 
results (31). The reason for the formation of abortive 
complexes could be either that there are too few Cre re- 
combinases at the loxP site or that there are one or more 
inactive Cre recombinases occupying the loxP site. 
However, an alternative interpretation that cannot be 
ruled out is that this 19.7% of the population of 
abortive complexes, results from inactive Cre molecules 
formed during preparation of the Cre enzyme. In spite 
of the observed formation of abortive complexes, 61 of 
76 molecules formed synaptic complexes rapidly from 
short-lived presynaptic complexes, and rate constants of 
(7.7 ± 0.1) x 10 _1 s" 1 and (4.2 ± 0.7) x 10" 1 s" 1 were 
obtained for the formation of wayward complexes and 
stable synaptic complexes, respectively. It has been 
reported that this process, which requires a collision 
between DNA partners, takes place relatively slowly 
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(31,46). However, rapid synaptic behavior has been 
observed in both the X Int and Cre recombinase systems 
(12,19,31). The similarity between these two tyrosine 
family recombinases reinforces the hypothesis that there 
are no intrinsic mechanistic properties that cause synapsis 
to be a step that occurs slowly in vitro. However, the dis- 
crepancy between the TPM results and in vivo observa- 
tions could be explained by the in vivo synapsis process 
being either intermolecular or involving a much larger 
separation between the recombining sites, making 
synapsis a slower process in vivo. 

Among the DNA molecules that formed synaptic 
complexes rapidly via a short-lived presynapse state, 25 
of 61 DNA molecules were subjected to a pathway 
ending in a wayward complex. These wayward complexes, 
in which the Cre-Cre interactions are relatively unstable, 
readily decompose at a rate constant of (2.7 ± 0.5) x 
10~ 2 s _1 . It has been shown that CreY324F is deficient 
in synapsis, showing a lower synapsis K d > 200 nM, and 
lacks in the strand cleavage (12). That only the wayward 
complex, decomposing at a rate constant of (3.2 ± 0.3) x 
10~~s _1 , was observed in CreY324F-/oxP system re- 
inforces the hypothesis that there are side pathways and 
futile intermediates in Cre-/oxP recombination process. 
Because of the limited resolution of TPM, it is not 
possible to differentiate whether the observed decompos- 
ition of these inactive synaptic complexes is actually due to 
the dissociation of Cre recombinase from one RBE of a 
single loxP site or is simply the decomposition of the 
wayward complex as a whole. 

Reversibility in strand cleavage/ligation is an expected 
property of the Cre-mediated site-specific recombination 
process because this reaction is isoenergetic, and no 
high-energy cofactors are required (4,25). If the reaction 
were irreversible, the synapse state duration analysis of 
DNA molecules containing parallel loxP sites would 
only reflect the decomposition of wayward complexes; in 
contrast, the synapse state duration analysis for DNA 
molecules containing inverse loxP sites will represent two 
reaction pathways, i.e. the decomposition of the wayward 
complex and the decomposition of the recombined 
synaptic complex. Based on kinetic analysis of the 
synapse state duration, similar rate constants of 
(2.0 ± 0.1) x 10 -3 s _I and (2.2 ± 0.2) x 10 -3 s" 1 were 
obtained for the rate-limiting step in the catalytic recom- 
bination process, i.e. the decomposition of the recombined 
synaptic complex, for DNA molecules containing parallel 
loxP sites and inverse loxP sites, respectively. All the 
experimental results support the hypothesis that strand 
cleavage/ligation is reversible in the Cre-mediated site- 
specific recombination process. 

This reaction pathway, strand cleavage/ligation, is 
abolished once the Tyr324 has been replaced with 
Phenylalanine, blocking the opportunity to entering the 
catalytic recombination process. An alternative explan- 
ation for the single reaction pathway observed in the 
decomposition of synaptic complexes is that CreY324F is 
simply defective in synapsis, and cleavage is not required 
for the formation of a stable synaptic complex. In order 
to verify the requirement of strand cleavage in the forma- 
tion of stable synaptic complexes, same experiments 



were performed in CreK201A-/oxP system. Based on 
kinetic analysis of the synapse state duration for 
the synapsis-capable CreK201A, even though it lacks in 
strand-cleavage ability, similar rate constant of 
(2.6 ± 0.2) x 10~ 3 s _1 , representing the decomposition of 
the synaptic complex, was obtained. It supported the state- 
ment that strand cleavage is not required for the formation 
of the stable synaptic complex in the Cre-loxP system. 

CONCLUSION 

Using a novel single-molecule method, TPM, the detailed 
reaction mechanisms and kinetics involved in the 
Cre-mediated site-specific recombination process have 
been revealed. The proposed reversible reaction mechan- 
ism and the presence of a thermally stable Holliday 
junction intermediate during the Cre-mediated site-specific 
recombination process are consistent with previous experi- 
mental results. Similar to the X Int system, wayward 
pathways and various reaction intermediates have been 
observed in the Cre-/oxP system too (19). Moreover, the 
detailed kinetic analysis performed in this study reinforces 
the importance of divalent cations and electrostatic inter- 
actions for stabilizing the close contact that maintains 
the Cre-Cre interface in the synapse state. These effects 
under different conditions were explored in detail using a 
single-molecule tool approach, which now offers a novel 
means of exploring reaction mechanisms other than con- 
ventional tools. 
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